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Glueballs at finite temperature from AdS/QCD 
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Inspired in the AdS/CFT correspondence, a variety of holographic phenomenological models have been proposed 
in the last years to describe non-perturbative aspects of strong interactions. These models are denominated as 
AdS/QCD. In this work we review the use of the AdS/QCD soft- wall model to investigate the spectrum of scalar 
glueballs at finite temperature. The scalar glueball states are identified as the poles of the retarded correlation 
function of the glueball operator. In the gauge/gravity duality, these poles are determined by the quasinormal 
spectrum of a massless scalar field propagating in the bulk geometry that consists on an AdSs black hole with a 
background dilaton field. We discuss some results for masses and decay widths of scalar glueballs in the plasma 
phase and analyse how these quantities evolve with temperature and momentum. 



1. Introduction 

The strong interactions are described by QCD 
which is a SU(3) gauge theory. At high energies, 
the coupling constant is small and the theory can 
be treated with perturbative methods. When one 
considers processes involving high energies one 
finds very good agreement between theoretical 
predictions and experimental results. However, 
many important aspects of strong interactions, 
like mass generation and quark confinement, are 
not yet completely understood because they in- 
volve the low-energy non-perturbative regime of 
QCD. So, it is important to look for alternative 
tools to deal with this non-perturbative regime. 
A traditional approach is to use lattice simula- 
tions that still do not furnish a complete de- 
scription of low-energy QCD. A more recent pro- 
posal is to construct phenomenological models for 
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strong interactions based on the AdS/CFT cor- 
respondence |ll2l3j . These models are presently 
known as AdS/QCD. The simplest model, known 
as hard wall, consists in using an AdS slice 
with a size related to an infrared cut-off in the 
gauge theory. This model was used to estimate 
hadron masses |4I5I6I7I8I9| . It presents confine- 
ment at zero temperature |10) and deconfine- 
ment at high temperatures [11| . Another in- 
teresting AdS/QCD model that leads to linear 
Regge trajectories is known as soft wall. In this 
model a non-uniform background dilaton field 
plays the role of an effective infrared cut-off. This 
model was also used to calculate hadronic masses 
|12|13|14] . The confinement/deconfinement ther- 
mal phase trasition in both models was studied 
in Refs. [15ll6ll7j . It was found that there is a 
Hawking-Page gravitational transition such that 
the space with a black hole is thermally favored 
at high temperatures, while below some critical 
(model-dependent) value, the thermally favored 
space is a pure AdS. 

In this paper we are going to review some recent 
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progress in using the soft- wall AdS/QCD model 
to compute the spectrum of scalar glueballs in 
the finite-temperature plasma phase of QCD [TH] 
(See also Ref. [T^). We consider the AdS black- 
hole spacetime for all temperature values. As 
commented above, for temperatures lower than a 
critical value the black hole is thermodynamically 
unstable and can be interpreted as a supercooled 
Yang-Mills plasma. 

The massless scalar field in the 5-d bulk is dual 
to scalar glueballs in the 4-d boundary theory and 
the black-hole quasinormal modes correspond to 
the poles of the retarded correlation function of 
the operator O — Tr(F^) that creates scalar glue- 
ball states. With a numerical analysis, we obtain 
the real and imaginary parts of the quasinormal 
mode frequencies as a function of the tempera- 
ture and the spatial momentum. The numeri- 
cal computations at high temperatures were per- 
formed using series solutions of the equation of 
motion. At low temperatures, we found problems 
with the convergence of the series so it was nec- 
essary to use an alternative approach. We have 
chosen the Breit-Wigner method which has been 
recently applied to compute black-hole quasinor- 
mal modes in Ref. We calculate the re- 
tarded Green's function of scalar glueballs and 
verify that their poles coincide with the quasi- 
normal mode frequencies. The imaginary part of 
the retarded Green's function gives the spectral 
function, which we also analyze in this article. 

Other recent results in scalar glueballs can be 
found in [21]. For a general review on glueballs 
see [221 . 



2. The soft- wall model at finite tempera- 
ture 

Our main objective in this paper is to review re- 
cent work [18] on the spectrum of scalar glueballs 
at finite temperature in the soft-wall AdS/QCD 
model. For this purpose, it is considered a dilaton 
field <i> = cz^ in a finite-temperature background 
that consists of an AdS black-hole spacetime with 
metric 

ds^ = e^-^ [-fdt^ + dx'dx' + /"^ dz"^] , (1) 



where / = 1 — (z/zh)^, A = —\n{z/L) and L de- 
notes the AdS curvature radius. The coordinate z 
is defined in the range < z < Zh and z = cor- 
responds to the boundary of the spacetime. The 
parameter z^ indicates the position of the event 
horizon which is related to the black-hole Hawk- 
ing temperature T by the relation Zh = I/ttT, 
where T also represents the temperature of the 
boundary field theory. 

The AdSs black-hole spacetime considered here 
is associated to three distinct phases, each one 
corresponding to a range of values of the tem- 
perature parameter T = nT/ ^Jc. The thermal 
competition between the black hole and thermal 
AdS space is such that there is a Hawking-Page 
transition between spaces at k, 2.38644. This 
transition is interpreted in the boundary field the- 
ory as a confinement/deconfinement phase tran- 
sition. The black hole spacetime is the dominant 
configuration for T > Tc, while the thermal AdS 
space is the dominant one for low temperature 
values. Hence, for T < Tc, the plasma will be in 
a metastable (unstable) phase corresponding to a 
positive (negative) sign of the black-hole specific 
heat Cbh- It was shown in Ref. [Hj that Cbh > 
for T 2 > 0.75 and Cbh < for f ^ < 0.75 . 

In the soft- wall model [T^], the dilaton back- 
ground field ^{z) — cz"^ interacts with the bulk 
fields through the replacement of the action in- 
tegrals: J d^x^—gC J d^x^—ge^'^C. In par- 
ticular, the action for the massless scalar field 4> 
is given by: 

^ = -7-? / d^'x^ge-" g''''dM(pdN4> , (2) 

where guN is the black-hole metric given by ([T]) 
and Kio is the ten-dimensional gravitational con- 
stant. The indices M, N run over 0, 1, 4 where 
x'^ (^ = 0, 3) are the 4d boundary coordinates 
and x'^ = z IS the extra radial coordinate. 



3. The Schrodinger equation 

We do not consider here the backreaction of the 
dilaton $ and of the scalar field on the metric. 
Then, from the action we obtain a linearized 
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equation of motion for 



After Fourier decomposition of (j) with respect to 
the coordinates x'^, Eq. ^ becomes 



(4) 



where A;^ = {-uj.qi), = and 

B = <S>~3A = cz^ + S \n{z/L). 

The foregoing equation takes the form of a one- 
dimensional Schrodinger equation when written 



in terms of the new variable ■(/' 
9^> + cj2i/, = y^, 



(5) 



where r* is the tortoise radial coordinate, defined 
in such a way that dz/dr^ — —f, and the effective 
potential V is given by 



V = I 



15 

T 



9z^ 



2cz^ 1 



(6) 



The general solution of equation ([5]) can be 
written as a linear combination of a normalizable 
solution ipi and a non-normalizable solution ^/;2. 
These wave functions have the following asymp- 
totic form near the boundary 2 = 0: 



i'2 = [1+0212^+0222'' + •••] 

-I- 6-01 ln(cz2), 

where the coefficients above are given by 

- g2 _ 2c 
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h = 
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(Cj2-g2)(^2_^2_4^) 

32 



(7) 
(8) 

(9) 
(10) 

(11) 
(12) 



The coefficient 022 is arbitrary and, in particular, 
we can choose 022 — 0. 



As it can be seen from Eq. © , the potential V 
vanishes when z ^ Zh so that Eq. ^ takes the 
form of a classical harmonic-oscillator equation 
with frequency uj. The solutions of such equation 
are proportional to exp(±ia;r*), being denoted by 
■0_ and -0+. The minus(plus)-sign solution can be 
interpreted as an ingoing (outgoing) plane wave 
at horizon z = Zh. The wave functions 'ip± can 
also be expressed in terms of the normalizable 
and non-normalizablc solutions: 



^± =-4(±)02+S(±)0i, 



(13) 



where A[±) and B(^±) are connection coefficients 
associated to Eq. (O, and determined as func- 
tions of UJ and q by imposing boundary conditions 
on the field solution. 

4. The glueball spectral function 

4.1. General expression 

The Minkowskian prescription of Rcf. 23J can 
be used to compute the glueball retarded Green's 
function G^'{uj,q) and, consequently, the corre- 
sponding spectral function TZ — —2lmG^{u!,q). 
First we write down the on-shell version of the 
scalar field action using the equation of mo- 
tion ©: 



S.,, 



d4xy^e-*5^^0a,0|„„, (14) 



So, we use the "bulk to boundary propagator" 
0fe(z) to decompose the scalar field 0(z, fc) as 
4>{z,k) = 0fc(z)0o(fc). In terms of 4>k{z), the ac- 
tion takes the form 



d^k 



0o(-fc).F(fc,z)0o(fc)L^„.(15) 



According to the Minkowskian prescription [51 
the retarded Green's function is then given by 



G"ik) = -2T{k,z 



„ „ lim ^ 



-99 



(16) 

'*0^(z)a,0,(z). 



where the propagator (j)k{z) satisfies an incoming- 
wave condition at horizon and the normalization 
condition lim.^o 4>k{z) = 1 at the boundary. 
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Figure 1. Left: Spectral function for g = and selected values of temperature. Right: Spectral function 
for ~ 0.05 and selected values of the momentum. The spectral functions are in units of N^/4tt'^. 



The function 4>k{z) can be expressed in terms 
of the wave functions defined in the last section 
as 



V'- 



A 



(-) 



3/2 i 

z ' e 



+ ^r^V-i 



A 



(-) 



(17) 

Then, using the relation t:^L^/2kIq = N^/Stt^ 
and the expansions ([7]) and ^ for tpi and 1/^2 , we 
obtain 

G^'W = ^ [(c + 2a2i)e-2 + 461n(ce2) 



2(5 + a2ic- 



+4 Re 



A 



i Ini- 



(-) 



(18) 



where e is an ultraviolet regulator. In the above 
equation, the e-dependent terms are ultraviolet 
divergent terms that can be removed using a holo- 
graphic renormalization procedure [24] . 

4.2. Numerical results 

To find the giueball spectral function we need 
to compute the imaginary part of the retarded 
Green's function G^{uj,q). From equation it 
follows that 



7^(.,,)^-,Im-^, 



(19) 



The ratio can be computed by a direct 

numerical integration of equation ([5]) following 
the steps described, for example, in Refs. [18l25j . 

In figure [U we show the dependence of the 
spectral function with the frequency uj — uj/^/c 
for some values of temperature T = nT/y^ and 
momentum q — q/^/c. The various peaks of 
Tl{u!,q) correspond to the poles of the retarded 
Green's function and are interpreted as glue- 
ball states in the gauge theory. According to the 
gauge/gravity duality, these poles are associated 
with the frequencies of the black-hole quasinor- 
mal modes. 

As it can be seen in the left panel of figure [l] 
the number of giueball states and their lifetimes 
(the inverse of the half-width) decrease with the 
temperature. In fact, there is a temperature value 
(T^ ~ 0.1) from which on there are no more peaks 
in the spectral function, characterizing a giueball 
melting. For a fixed value of temperature T, the 
position and the width of the peaks increase with 
the momentum q, as shown in the right panel of 
figure [H 

5. From quasinormal modes to glueballs 

Two different methods have been used to com- 
pute the scalar-field QNM frequencies of the AdS 
black-hole spacetime ([T]). The first one is a series 
expansion method [26j , which reduces the prob- 
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Figure 2. Numerical results for the square of the real and imaginary parts of the QN frequencies, oj'^/c 
and wf/c, for the first five quasinormal modes n = 0, 1, 4, with q = 0. 



lem to finding roots of a polynomial, and it is 
suitable to find QNM frequencies in the high and 
intermediate temperature regimes. The second 
method is based on the computation of Breit- 
Wigner ressonances [20], and it is efficient to in- 
vestigate the spectrum for very low temperatures. 
The approaches are complementary and present 
a very good concordance in a wide range of pa- 
rameter space P^ . 

We show in Fig. [2] the real and imaginary parts 
of the QNM frequencies as function of the tem- 
perature for the first five modes n = 0, 1, 4 
with 5 = 0. In the limit T ^ 0, the quasinor- 
mal modes tend to the spectrum of glueballs at 
zero temperature: =4c(n-|-2). As the tem- 

perature increases the imaginary part of the fre- 
quency uji also increases, while the real part uiji 
decreases with T. It means that the mass and 
the decaying^ time of the glueball excitations de- 
crease with T for very low temperature values. In 
the intermediate- and high-temperature regimes, 
the effect of the dilaton background field is small, 
and the QNMs are essentially the same as those of 
the black hole solution with $ = j27|28|29|30] . 
In particular, after reaching some minimum value 
the real part ujr increases with temperature tend- 
ing asymptotically to a linear dependence. In 
fact, the fractional differences between the cases 
with and without dilaton reach a value of the or- 
der of 0.15 for T = I, and tend asymptotically to 



zero in the limit of high temperatures [18j . 



6. Final comments 

In this article we review recent progress [18] on 
the study of the spectrum of scalar glueballs at fi- 
nite temperature. The AdS/QCD soft-wall model 
and a black-hole background spacctimc have been 
used to investigate thermal effects on glueball ex- 
citations of the plasma phase. As it was shown 
in Ref. [T^ the black-hole space are associated 
to an unstable, metastable or stable phase of the 
dual theory according with the value of the tem- 
perature parameter T = T / ^/c. It was noticed 
the existence of peaks in the spectral function, 
corresponding to glueball states, only for temper- 
atures lower than = 0.1, a regime for which the 
black hole is in an unstable phase. So according 
to our results, we do not expect to find glueballs 
in the hot plasma phase (T^ > 2.39) of the QCD. 
Similar results were recently obtained for scalar 
[19] and vector [31] mesons in the soft-wall model. 
It would be interesting to extend these computa- 
tions to other fields and other AdS/QCD models, 
like those of Refs. 32 33 34J. Such analysis could 
reveal what aspects of our results are specific of 
the holographic model considered here and what 
aspects represent general characteristics of glue- 
ball excitations in the dual plasma. 
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